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Abstract 

The idea of isotropic resonant magnetism in the visible range of frequencies known from prece- 
dent publications is developed having in mind achievements of the modern chemistry. Plasmonic 
colloidal nanoparticles covering a silica core form a cluster with resonant and isotropic magnetic 
response. Two approximate models giving the qualitative mutual agreement are used to evaluate 
the magnetic polarizability of the cluster. It is shown that the electrostatic interaction of nanocol- 
loids decreases the resonant frequency of an individual complex magnetic scatterer (nanocluster) 
compared to the previously studied variant of a planar circular nanocluster with same size. This 
means the reduction of the optical size of nanoclusters that presumably allows one to avoid strong 
spatial dispersion within the frequency range of the negative permeability. 

PACS numbers: 78.20.Ci, 42.70.Qs, 42.25.Gy, 73.20.Mf, 78.67.Bf 
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I. INTRODUCTION 



The interest in artificial resonant magnetism in the optical range which was recently 
inspired by the development of metamaterial science (see e.g. in [l[ Q]) is expected to 
grow after successful demonstration of the subwavelength optical imaging obtained in the 
far zone of objects with the use of so-called hyperlenses 3|, |4|. Hyperlenses transform 
evanescent waves into propagating ones and create the spatially magnified optical images 
keeping subwavelength details of the objects. 

Known hyperlenses operate with TM-polarized waves. This polarization is related to 
the negative permittivity of metals (which are used by known hyperlenses) in the visible 
frequency range. Similar operation with the TE-polarized light would require materials 
with negative permeability. 

Materials with negative fi in the visible range do not exist in nature. Recently, many 
attempts to engineer such materials were reported. They include, for example, lattices of 
paired pla smonic nanowires, nanoplates, nanocones, fishnets, and plasmonic split rings j^, E3, 



3,0,3, 
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121 ]. These structures are all geometrically strongly anisotropic. The strongest 



magnetic response corresponds to one specific direction of the propagation of electromagnetic 
wave. Notice, that in cited papers the magnetic response of these structures was reported 
only for this single direction. The possibility to use strongly anisotropic artificial magnetic 
media in superlenses and hyperlenses has not yet been studied. 

The first work devoted to the possibility of the isotropic optical magnetism that cannot 
be achieved in the mentioned anisotropic structures was, to our knowledge, [l^J]. The idea 
of the present paper is to link this resonant magnetism to the existing technology and to 
find the design which ensures the real isotropy of the resonant permeability, i.e. absence 
of strong spatial dispersion. We use the term strong spatial dispersion for the case when 
the dispersion diagram of the lattice is quantitatively different from the dispersion in any 
continuous media. Weak spatial dispersion is a special term applied to designate effects of 
artificial magnetism and bianisotropy in composite media [l^. 



In 



14| one suggested to prepare optical magnetic media using plasmonic (colloids of 



noble metals) nanospheres arranged in nanorings. In high-frequency magnetic field H 
nanocolloids are polarized in the azimuthal direction with respect to this field. To obtain 
an isotropic resonant magnetic medium one should fabricate a dense array of nanorings 



lying in 3 mutually orthogonal planes. The dense isotropic package means that the array 
of nanorings must be a fee or a bee lattice, that allows us to consider the structure as a 
phobic crystal. The theoretical prediction of fl shows that in this case the _e of 
the effective permeability /x e ff can be strong enough so that the real part of fi c s takes negative 
values. However, from [14| it has not yet been clear how to fabricate these structures and 
their isotropy was not proved. In fact, an isotropic artificial magnetism is not possible for 
photonic crystals (lattices with strong spatial dispersion). 

Formally, the resonant magnetic permeability as well as the permittivity can be attributed 



to photonic crystals as well (see e 



. g . Q,0 : 



17]). However, at frequencies where a cubic 



lattice becomes a photonic crystal the refraction index and the wave impedance both become 
strongly dependent on the propagation direction [22[]. This implies that effective material 
parameters e e s or/and /i fr of the crystal are very different for a wave propagating along 
one of the lattice axes and for a wave tilted with respect to this axis. This evidently means 
the absence of the isotropy. Of course, the random mixture of particles is more optically 
isotropic than the lattice with the same concentration of particles. However, in the present 
case the requirement of the negative magnetic response of the composite medium obviously 
implies the very dense packaging of magnetic scatterers. In the present case we cannot 
suppress the spatial dispersion using the randomization of the array. 

The comparison of the approximate dispersion diagram for the homogenized lattice with 
the exact dispersion diagram was done in some papers, e.g. in 2CJ? ]. In 20| the fee (face- 
centered cubic) and bec (body-centered cubic) lattices of metal-covered dielectric spheres 
were studied. It was found that in fee and bec lattices the spatial dispersion effects arise 
within the frequency band of the plasmonic resonance. The spatial dispersion in this region 
leads to closing the resonance band-gap with simultaneous strong anisotropy of the lattice 
refraction index. These effects correspond to the frequencies at which the distance between 
adjacent particles centers exceeds one quarter of the wavelength in the matrix. 

Unfortunately, in the known literature there are no general results for the frequency 
intervals where the effects of strong spatial dispersion arise in lattices. Of course, it is clear 
that the frequency range corresponding to D/X > 0.5 definitely implies the strong spatial 
dispersion since A = 2D is the Bragg resonance condition. In lattices of strongly scattering 
inclusions the condition of the absence of strong spatial dispersion becomes more restrictive 
since the polarization of inclusions shortens the effective wavelength A e fr compared to A. 



It is known that the artificial resonant magnetism in arrays of small complex-shape metal 
particles at microwaves practically corresponds to relation D p /X r > 0.1, where D p is the 
typical size of the particle 13j, |21[ and A r is the wavelength of the magnetic resonance of a 
particle in host medium. Notice, that the ratio D p /X ~ 0.1 . . . 0.15 is practically achievable, 
e.g. for well-known double split-ring resonators (SRRs). SRRs of this optical size can be 
considered as small enough to avoid the strong spatial dispersion in a compact lattice. Then 
we can consider the low-frequency stop-band of this lattice associated with the resonance of 
SRRs as the effect of the negative permeability. Really magneto-electric media are opaque 



when their material parameters have different signs / 



However, the metamaterial suggested in [14j corresponds to a quite large optical size of 
a magnetic scatterer at its resonance, where the condition D p /X r < 0.15 is far from being 
fulfilled. Really, consider the results for /i e g- illustrated by Fig. 2 (a,b) of [14J . The negative 
permeability was theoretically achieved at 630 or 680 THz with nanorings of the following 
geometry. Six or four metal nanospheres of radius a p = 16 nm are centered at the circle of 
radius Ro = 40 or i?o — 38 nm. This means that the size D p of the effective magnetic particle 
(i.e. the distance between the edges of opposite metal nanospheres) is 112 or 108 nm. The 
wavelength in the host medium with background permittivity = 2.2 at 630 THz is equal 
to X r = 323 nm. Then the size of the nanoring from 6 spheres is D p pa 0.35A r i.e. the optical 
size is rather large. The packaging of these nanorings in the host medium was assumed in [ijj] 
to be very dense, namely corresponding to their concentration N^r = (108)~ 3 nm -3 . This 
concentration is physically achievable for a bcc or a fee lattice of nanorings. However, even 
so dense packaging corresponds to the distance between the centers of adjacent nanorings 
D pa 0.4A r . Having in mind results of 2^] we believe that in this case the strong spatial 
dispersion arises and the isotropic magnetism is impossible. 

The clear way to decrease spatial dispersion in a lattice of magnetic scatterers i.e. to 
suppress the angular dependence of /x e ff is to decrease the size of a scatterer keeping its 
magnetic polarizability. In this way we will also decrease the optical distance between the 
centers of adjacent magnetic particles. Below we modify the design of magnetic particles 



suggested in 



14j in this way. We obtain the ratio D p /X r ~ 0.2 . . . 0.25 which is presumably 



below the range of the strong spatial dispersion. Simultaneously, we link our design to 
existing nanotechnologies. 
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II. CORE-SHELL MAGNETIC CLUSTERS 



The design of the optical magnetic scatterer suggested in this paper is based on the 
possibility to fabricate arrays of silver or gold colloidal nanoclusters located on spherical 
dielectric cores in a liquid host medium. The fabrication technology uses self-assembled 
adhesion of metal nanospheres on the silica core. The idea of such self-assembling aggregates 
was probably first realized in 



23]. However, this process led to the complete covering of 
the silica core by contacting gold nanoparticles. The process chemically stimulated by 
some acids which allowed one to control its speed and to obtain the needed number of 



nanocolloids per one core was described in 



241 ]. Later, nanoclusters comprising few tens of 



colloidal particles on a silica core of diameter of the order 100 nm were reported in works 
25)]. For obtaining the maximal electromagnetic response for given sizes (of the core and 
of the nanocolloids) their number N tot per one nanocluster should be as large as possible. 
However, it is difficult in this way to prevent the touching of nanocolloids which is destructive 
for their plasmonic resonance that will degrade due to the tunneling of electrons between 
nanocolloids. Practically, minimal separation d between nanocolloids should be comparable 
with the radius of colloidal particles a p . Below we respect the condition d ~ q „/2 that 
restricts N tot . This can be practically obtained using the technology described in [26( and 
2?]]. In these works the silica core of diameter 100 nm (the technology allows in principle 
to reduce this size to 20-30 nm) was covered by mutually touching polystyrene nanospheres. 
Instead of simple polystyrene nanospheres one can use core-shell particles with plasmonic 
nanocolloids as cores. This method will allow us to reliably control the separation d between 
colloidal particles which is simply equal to the double thickness of the polystyrene shell. 

Such a nanocluster is shown in Fig. [1] (a) (a general view) and (b) (a cross section). It 
should have the isotropic resonant response to the electromagnetic field of light. At one 
frequency the induced electric dipole dominates, and the polarization of colloidal particles is 
parallel to the electric field of light. At another frequency the induced magnetic dipole dom- 
inates, and the polarization of colloidal particles is azimuthal with respect to the magnetic 
field. In other words, the applied magnetic field form the effective nanorings of colloidal 
spheres glued to the silica core. Both electric and magnetic resonant frequencies origin 
from the plasmonic resonance of an individual colloidal nanosphere. In this paper we study 
only the magnetic resonance of such a nanocluster. Therefore we call it a magnetic nan- 





(a) (b) 

FIG. 1: A magnetic nanocluster (a bulk isotropic scatterer with the magnetic resonance): (a) - a 
general view of a cluster of metal nanocolloids with dielectric core, (b) - the resonant azimuthal 
electric polarization is excited in the discrete "shell" formed by plasmonic nanocolloids when 
the applied field is the optical magnetic field directed along z. 



ocluster (MNC). Bulk arrays of MNC can be created not only in liquid but also in porous 
matrices where the separation between the pores should be small enough to obtain the high 
concentration of MNC and a strong magnetic response of the medium. 

Below we will see that the frequency of the magnetic resonance of MNC is strongly 
reduced compared to the frequency Q p of the plasmonic resonance of a single nanocolloid. 
This effect results from the strong electrostatic coupling between the nanocolloids in the 



14[ , and it leads 



MNC. This electrostatic coupling exists also in a single nanoring studied in 
to a certain reduction of the magnetic resonance frequency compared to Q p . However, in 
the present geometry this reduction is much stronger since the most important electrostatic 
interaction is that between colloids having parallel electro-dipole moments. 

The total number N tot of colloidal nanospheres per one MNC can be controlled by the 
self-assembling technological process. Also, it can be adjusted by properly chosen initial 
concentrations of nanocolloids and silica spheres in the host liquid. N tot can be expressed 
through the core radius a, the colloid radius a p and the separation d between metal spheres 



(a) 



(b) 



FIG. 2: Two models of the magnetic polarization of MNC in the applied magnetic field directed 
along z: (a) - nanocolloids form regular rings around the z— axis, the colloidal particles are equidis- 
tant in the rings of radiuses Rj (induced dipole moments are shown by bold arrows) and the rings 
are equidistant; (b) - nanocolloids form around the core (shown as a dashed circle) an effective 
continuous shell with azimuthal bulk polarization Pa,. The shell of thickness 2a p + d can be char- 
acterized by resonant isotropic permittivity e^. 



Here [A] denotes the integer part of the number A. In this formula the curvature of the 
portion A of the spherical surface with radius Rq = a + a p + d/2 which is cut of the sphere of 
radius Ro (this sphere centers the nanocolloids and their dipole moments are located on it) 
by one core-shell particle of radius a p + d/2 is neglected. In other words, in our calculations 
the surface A of the effective sphere Rq per one colloidal particle is assumed to be a planar 
square with side 2a p + d. 

III. TWO MODELS OF MNC 

Two approximate models for obtaining the magnetic polarizability a mm of a MNC are 
illustrated by Fig. [2] (a) and (b). The agreement of their results can be hopefully considered 
as a validation of these results. These models describe in two different ways the electromag- 
netic interaction of colloidal particles. In the first model the nanocolloids on the dielectric 



(see Fig. HD: 




(1) 
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core are assumed to form regular rings around the z— axis. The rings of radiuses Rj are dis- 
tanced by a p /2 from one another. The distance between the nanoparticles in any ring is also 
equal to d — a p /2. The presence of the polystyrene shell of nanocolloids does not influence 
to the result, since the permittivity of the shell is equal to that of the matrix (e n = 2.2). 

In the second model the discrete structure of core-shell plasmonic nanoparticles on the 
central silica core is treated as a nanolayer of effective medium of thickness 5 = 2a p + d 
covering the silica core. The azimuthal electro-dipole polarization and the isotropic 
permittivity El of this layer can be easily calculated and we find a mm in the closed form 
calculating the magnetic moment of this homogenized shell. 

The presence of the central silica core also does not influence to the result due to the 
absence of its electric polarization by the applied magnetic field and smallness of its magnetic 
moment which is not resonant. It was checked analytically that the interaction of this 
small magnetic moment with electro-dipole moments of colloidal particles is negligible. The 
permittivity of the silica core e s = 4 > e n = 2.2 influences only weakly to the response of 
colloidal spheres. The polarizability a of a single colloidal particle is assumed to be the same 
as if nanocolloids were located in the uniform medium with permittivity e n '- 



a 



4irals e h - 



(2) 



"e m + 2e h J 6n n 
Here e m is the permittivity of the metal which is taken in the same form as in jl4j ]: 

F F (27T/ p ) 2 pKfUU 

<^m — fc i o r * 5 • \o ) 

U> 2 US 1 

For silver colloids we have the plasma linear frequency f p = 2175 THz, the damping fre- 



quency fa = 4.35 THz and t 



re ultraviolet permittivity = 5 



Q. 



As it was explained in [14J the magnetic moment of nanorings is excited by the magnetic 
field of the light wave. To calculate the magnetic polarizability of a MNC we assume that 
it is excited by magnetic field polarized along z [see Fig. [TJ (b)] that in the model of the 
magnetic excitation assume to be uniform across the circle of radius Rq and varying with 
time as exp(— tut)) with amplitude H Q . From Maxwell's equations this applied magnetic 
field H a is associated with the azimuthal electric field E a vanishing at the center of MNC: 

H a = z H , E a = <f (ir}H kR), (4) 

where rj = a/ Jio/so^h 1S the wave impedance of the host medium of permittivity s n - Unit 
vectors zq and </>o correspond to the cylindrical coordinates (R, <f), z) shown in Fig. [1] (b). In 



fact, the applied electric field given by (J4]) is associated with slightly non-uniform magnetic 
field H a = z H (l — k 2 R 2 /2), but this small correction to the applied magnetic field matters 
nothing for our model. We calculate the magnetic moment of MNC induced by H a whose 
amplitude at the center of MNC is H through the electric polarization of nanocolloids by 
the electric field associated with H a . 

The 0-oriented dipole moments pj of colloidal particles of any j-th ring are identical 
in polar coordinates. The z-directed magnetic moment rnj of the j-th ring can be found 
using formula (8) of [l4]. The total magnetic polarizability of MNC is then obtained as 



a m m = Y^ a mm, where a™L = mj/H . 

Notice, that the obvious (not sufficient) condition of the absence of spatial dispersion in a 
bulk array of MNC is the independence of the magnetic polarizability a mm of an individual 
MNC on the distribution of the applied magnetic field around the center of MNC, if the 
condition of the magnetic excitation is respected: H a (R = 0) = z H , E a (R = 0) = 0. This 
condition can be satisfied for example by a standing plane wave, i.e. a pair of plane waves 
travelling along the a;— axis: 

H a = z ^-(e-i kx + e***), E a = y^{e~ lkx - e jkx ). 

The comparison of a mm calculated for this azimuth-dependent excitation of MNC and that 
calculated for excitation (Tj0) would be a useful exercise, and we plan it for the next paper. 
It could indicate the frequency bounds of the spatial dispersion related to the finite size of 
MNC. However, in the present paper we restrict the analysis by a more simple excitation 
case @. 

In both models the reduction of the resonance frequency of MNC compared to the plas- 
monic resonance Q p of the individual colloidal particle in the host medium with Eh = 2.2 is 
determined by the electromagnetic coupling of colloids. In fact, this reduction will be even 
stronger than it is predicted below. This is because the presence of the silica core makes the 
effective host medium permittivity slightly larger than e^. We neglect this small effect for 
simplicity of the model, however, its role is positive. 



A. The model of regular rings 



Consider a dipole p = p(j)4>o shown in Fig. [T] (b) that belongs to the j-th effective ring 
of MNC. Assume that N tot ^> 1 (practically the model is applicable when N tot > 15 . . . 20). 
Then it is easy to show using auxiliary spherical coordinates and formula ([T]) for N tot that the 
radius of the j-th ring is approximately equal to Rj ~ Rq sin(2j ^/7r/N tot ). The z— coordinate 



of the j-th ring is z(j) ps R cos(2j^n/N tot ). The number of dipoles Nj in the ring is equal 
Nj 2TcRj/(2a p + d). The number of such rings in MNC is equal N r ps ^TiN tot /2. The 
cylindrical coordinates of the dipole p in the j-th ring are z = zj, R = Rj, = 2nq/Nj where 
q = . . . (Nj — 1) is the integer number determining the position of the dipole within the 
ring. 

The field produced by the dipole p at the point with coordinates (z f , R', 0') can be found 
from the standard formula: 



E(R, R') 



e ~jkr / e ikr ike ikr 

k 2 (r x p) x r — h (3rp ■ r — pr 2 ) I — — 



^3 \ i^*4 



(5) 



Ane e h 

Here r = ^/R 2 + R' 2 - 2RR' cos(0 - 0') + (z — z') 2 is the distance between radiating p(j) 
and receiving p(ri) dipoles of MNC. Due to the azimuthal symmetry of the problem only the 
0— th component of the field produced by all rings of dipoles p(j) is nonzero at the center 
of the receiving dipole p(n), i.e. at point (z', R', 0') . In other words the local electric field 
acting on the dipole p(n) is orthogonal to the vector R! shown in Fig. [T] (b). Respectively, 
only the azimuthal component of the vector E(R, R'), should be taken into account. 

From (jSj) the scalar interaction coefficient of dipoles p(j) and p(n) can be easily derived. 
It is defined as the 0-th component of the field produced by the unit azimuth-oriented 
dipole located at point (z = Zj,R — Rj, = 2irq/Nj) and calculated at point ( R! = 

i? n ,0' = 2vr S /A„): 

n ikr 



Q n l = 



(kr) 2 (R' 2 + (z- z')R' - RR' + cos(0 - <p')(R 2 + (z - z')R) - RR 1 ) 
3 J R J R'sin 2 (0-0')(l - ikr) -cos(0-0')(l - ikr)r 2 ] . (6) 



The dipole moment of the receiving dipole is equal to p — aE loc . Here the local field is the 
sum of the external electric field (0J E a (R') = (iH r]kR f /2) and all dipole fields: 



E loc = E a (R')+p(j)J2Q 



>> j 

qs ■ 



in 



This way we obtain the system of equations for dipole moments of colloidal nanospheres of 
any ring: 

9=1 j^n 9=0 

Here the term with Qq™ describing the interaction of the dipole p(n) located at (z 1 = z n , R' = 
R n , (J)' = 0) with other dipoles of the same n-th ring is shared out. The expressions for 
coefficients Qq™ entering this term and given by (J6j) with z = z' are described by (EH). There 
is no coefficient Qqq in (j7|) since the summation starts from q = 1 (see also in 141]). 

Solving the system (171) we find all dipole moments p(n). The magnetic moment of j-th 
ring is calculated as in [14 1: 

m . = -^(j)a^ (g) 

Then we obtain the magnetic polarizability of the MNC as the sum of m~ letting H = 1: 



j 

N r 

2 ~ 



limn, — — £ ;PU) N jRj- 
3=1 



The relative effective permeability of the composite medium is given by 14j : 



1 

^ = 1 + Api -i ~ - (9) 

ls MNC U mm 3 

Here Nmnc is the volume concentration of MNC that can be expressed though the effective 
volume per one magnetic scatterer N MNC = 1/Vo = 1/D 3 . In simple cubic lattices D > 
Dp = 2(a + 2a p ), and D is the unit cell size. However a more compact arrangement of MNC 
is also possible when D < D p and N MNC > V-^p- 

The inverse polarizability of a nanocolloid in (jTJ) corresponds to formula (j2J) and contains 
the term (— ik 3 /67reoEh) that describes the radiation damping. The radiation damping of the 
magnetic dipole with magnetic polarizability a mm should be described by the term (— z/c 3 /67r) 
[14J]. It is known that the radiation damping is cancelled out in regular 3D arrays and for 
regular lattices we should have used instead of (jUJ) the relation [l4|: 

^ = 1 + AT' 1 fa-^ + ^-l ' (10) 

^ MNC \ a mm ' 1 Q n ) 3 

However, the dissipative losses due to the plasmonic resonances of metal nanospheres 
strongly dominate over the contribution of radiation losses into the imaginary part of the 
permeability, i.e. the difference in results of (jUJ) and of ffTUl) is negligible. 

1 1 



B. The model of an effective shell 



Replacing the discrete plasmonic structure by the continuous shell we introduce the bulk 
polarization P^ that can be expressed through the averaged 4>— polarized electric field E &y 
in this shell and its unknown permittivity e^: 

P<p = e (e L - e h )E av (11) 

The field E av is related with the local field acting on any colloidal nanosphere by the Clausius- 
Mossotti relation: 

E av = E\ oc — — — . (12) 

Here V\ = 2a p (2a p + d) 2 is the volume per one colloidal nanosphere and p = P^Vi is the 
dipole moment of the reference nanosphere. Using the formula p = aE\ oc together with ifTTl) 
and ( fl2l) we come to the Lorentz-Lorenz formula: 

£ L = £ hi 1 + aeoel l j ■ ( 13 ) 

The definition of the magnetic moment of any volume V comprising polarization currents 
j reads as: 

1 

m = - / j x rdV. 
2 J J 

v 

It can be rewritten for the MNC in terms of the bulk polarization Pf. 



m = j P^RdV. (14) 

v L 

Here Vl = 4:nRl(2a p ) = 87ra p (a + a p ) 2 is the volume of the spherical layer with central radius 
Rq and the thickness 2a p . The integration of the bulk polarization across this layer can be 
replaced by simple product P^(2a p ). The polar radius R that enters (JT3J) can be expressed 
in spherical coordinates as R{9) = Rq sin 9. After substitution of ( fTTj) into ( fl4l) we obtain: 

2-7T 7T 

m = -icu£ a p (e L - e h )Rl J d<f) J d6 sin 9E av (9)R(9). (15) 

o o 

The averaged electric field of azimuthal polarization at the circle of radius R(9) is re- 
lated to the magnetic field H at the center of the MNC as E av (9) = iu){i R(9)H /2 = 



19 



(iHof]kRo sin 9/2), and after this substitution into (fl5|) and trivial integration we come (let- 
ting H — 1) to the following formula: 

a mm = m= —(k R ) 2 V (e L - e h ), (16) 

where it is denoted Vq = R^a p and k = k/ y/e~h = uo^/Eq^o is the free space wave number. 

Substituting ( f!3l) into (fl6l) we obtain the final closed-form formula for the magnetic 
polarizability of individual MNC: 

a mm = 47r(fco# )Vo - 3£o£ ^ (17) 
The permeability can be then found using ([91) or ( flOl) . 



IV. RESULTS AND DISCUSSION 

First, the results obtained in [jjj] were reproduced using the first model (the second model 
implies very large N tot and cannot be applied to the case of a single ring in MNC). Namely, 
the system ([7]) was solved for N r = 1, when it was taken also N tot = 4, Rq = 38 nm and 
a p = 16 nm. This geometry corresponds to the total size of the magnetic cluster D p = 108 
nm and to the spherical core radius (recall, that in the present theory the difference of the 
core permittivity from that of the host medium is not taken into account) a = Rq — a p = 22 
nm. The concentration Nmng of effective magnetic scatterers was assumed in [yj] to be 
equal Nmng — 95~ 3 nm -3 (almost touching MNC in a fee or a bec lattice), the metal of 
colloidal particles was silver with permittivity (J3]). 

The result for the effective permeability is presented in Fig. |3] (a). It reproduces Fig. 2 
(b) of [ijj] with high accuracy. This agreement can be considered as a check. 

If we keep the same separation of colloidal particles d = 10 nm as in the previous example 
the maximal possible number of colloidal particles for the spherical core of radius that 
corresponds to this example is equal to N tot = 16. It corresponds to N r = 3. Keeping the 
same concentration Nmnc — 95~ 3 nm~ 3 of same MNC we obtained with this geometry the 
result shown in Fig. [3](b). This result was also obtained using the first model (for N tot = 16 
the second model is inadequate). The reduction of the magnetic resonance frequency due 
to the presence of two additional nanorings compared to Fig. |3] (a) is not dramatic though 
visible. 
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FIG. 3: Effective permeability of the array of MNC with sizes D p = 108 nm, a = 22 nm, a p = 16 
nm (Ro = 38 nm) hosted in the matrix with Eh = 2.2. |(a)| One ring of nanocolloids in every MNC, 
the concentration of MNC is Nmnc = 95 -3 nm~ 3 . |(b)| Three nanorings in every MNC, other 
parameters are the same. Real and imaginary parts of the permeability are shown by solid and 
dashed lines, respectively. 



More promising results were obtained keeping the same total size of MNC D p = 108 nm 
for a p = 13 nm, d = 6 nm, that implies a = 22 nm. This geometry corresponds to one 
effective ring with 9 colloids, two rings with 5 colloids and two rings with 3 colloids. In this 
case the number of colloidal particles is large enough and the qualitative agreement of two 
models was obtained. The concentration of MNC in Fig. H] is the same as in the previous 
example: Nmnc — 95~ 3 nm' 3 . The two models show approximately the same resonant 
frequency and the same magnitude of the resonance. The second model ignores, of course, 
the higher magnetic resonances that should arise in interacting nanorings and are revealed 
if we use the fist model. However, these higher resonances are weak and have no practical 
meaning. 

The reduction of the resonant frequency in Fig. |4] compared to Fig. [3] means that the 
resonant size of the MNC reduces from 0.35A r as in the previous example to 0.25A in the 
present one. Here A r is the resonant wavelength in the host medium. 

High magnitude of the Lorentz resonance in Fig.HJmeans that we can also reduce the con- 
centration of MNC making the array of MNC easier for manufacturing. The results depicted 
in Fig. [5] (a) and (b) are obtained using the first model for Nmnc — 110~ 3 nm -3 instead of 
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FIG. 4: Effective permeability of the array of MNC with sizes D p 
nm, d = 6 nm hosted in the matrix with Eh = 2.2, the concentration of MNC is Nmnc = 95 -3 
~ 3 . (a) : First model of the interaction of nanocolloids. |(b)| Second model of the interaction of 



nm 



nanocolloids. 
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FIG. 5: Effective permeability of the array of MNC with sizes D p = 108 nm, a = 22 nm hosted 



in the matrix with = 2.2, the concentration of MNC is Nmnc = HO nm 3 . (a} 
d = 6 nm. |(b)[ a p = 12.5 nm, d = 6 nm. 



a p = 13 nm, 



Nmnc = 95~ 3 nm -3 as in the previous case. The difference between two plots in Fig. [5] (a) 
and (b) is determined by a 1 nm difference in the radiuses of plasmonic nanospheres, and 
demonstrates how sensible is the magnetic response of MNC to the deviation of parameters. 
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V. CONCLUSIONS 



In the present paper we developed the idea of the resonant optical magnetism in its 
isotropic variant modifying the known design of optical magnetic scatterers suggested in jjjl . 
The array of isotropic magnetic nanoclusters proposed in the present paper can be obtained 
in a liquid or porous matrix using the existing nanotechnologies. The metal nanocolloids 
are located in the core-shell particles attached to the silica core. Their electrostatic mutual 
coupling helps to reduce the resonant frequency for the given size of a nanocluster compared 
to the frequency of the plasmonic resonance of a single nanocolloid. This reduction is the 
essential result of the present paper since it allows one to decrease the spatial dispersion of 
the dense array (lattice) of nanoclusters. In other words, the suggested geometry is helpful 
for the isotropy of the effective resonant permeability of the array. The further optimization 
of the structure, further validation of the suggested model by full-wave numerical simulations 
are planned in the next paper. 
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Abstract 

The idea of isotropic resonant magnetism in the visible range of frequencies known from prece- 
dent publications is developed having in mind achievements of the modern chemistry. Plasmonic 
colloidal nanoparticles covering a silica core form a cluster with resonant and isotropic magnetic 
response. Two approximate models giving the qualitative mutual agreement are used to evaluate 
the magnetic polarizability of the cluster. It is shown that the electrostatic interaction of nanocol- 
loids decreases the resonant frequency of an individual complex magnetic scatterer (nanocluster) 
compared to the previously studied variant of a planar circular nanocluster with same size. This 
means the reduction of the optical size of nanoclusters that presumably allows one to avoid strong 
spatial dispersion within the frequency range of the negative permeability. 

PACS numbers: 



1 



I. INTRODUCTION 



The interest in artificial resonant magnetism in the optical range which was recently 
inspired by the development of metamaterial science (see e.g. in [l[ Q]) is expected to 
grow after successful demonstration of the subwavelength optical imaging obtained in the 
far zone of objects with the use of so-called hyperlenses 3|, |4|. Hyperlenses transform 
evanescent waves into propagating ones and create the spatially magnified optical images 
keeping subwavelength details of the objects. 

Known hyperlenses operate with TM-polarized waves. This polarization is related to 
the negative permittivity of metals (which are used by known hyperlenses) in the visible 
frequency range. Similar operation with the TE-polarized light would require materials 
with negative permeability. 

Materials with negative fi in the visible range do not exist in nature. Recently, many 
attempts to engineer such materials were reported. They include, for example, lattices of 
paired pla smonic nanowires, nanoplates, nanocones, fishnets, and plasmonic split rings j^, E3, 



3,0,3, 



10 



11 



121 ]. These structures are all geometrically strongly anisotropic. The strongest 



magnetic response corresponds to one specific direction of the propagation of electromagnetic 
wave. Notice, that in cited papers the magnetic response of these structures was reported 
only for this single direction. The possibility to use strongly anisotropic artificial magnetic 
media in superlenses and hyperlenses has not yet been studied. 

The first work devoted to the possibility of the isotropic optical magnetism that cannot 
be achieved in the mentioned anisotropic structures was, to our knowledge, [l^J]. The idea 
of the present paper is to link this resonant magnetism to the existing technology and to 
find the design which ensures the real isotropy of the resonant permeability, i.e. absence 
of strong spatial dispersion. We use the term strong spatial dispersion for the case when 
the dispersion diagram of the lattice is quantitatively different from the dispersion in any 
continuous media. Weak spatial dispersion is a special term applied to designate effects of 
artificial magnetism and bianisotropy in composite media [l^. 



In 



14| one suggested to prepare optical magnetic media using plasmonic (colloids of 



noble metals) nanospheres arranged in nanorings. In high-frequency magnetic field H 
nanocolloids are polarized in the azimuthal direction with respect to this field. To obtain 
an isotropic resonant magnetic medium one should fabricate a dense array of nanorings 



lying in 3 mutually orthogonal planes. The dense isotropic package means that the array 
of nanorings must be a fee or a bee lattice, that allows us to consider the structure as a 
phobic crystal. The theoretical prediction of fl shows that in this case the _e of 
the effective permeability /x e ff can be strong enough so that the real part of fi c s takes negative 
values. However, from [14| it has not yet been clear how to fabricate these structures and 
their isotropy was not proved. In fact, an isotropic artificial magnetism is not possible for 
photonic crystals (lattices with strong spatial dispersion). 

Formally, the resonant magnetic permeability as well as the permittivity can be attributed 



to photonic crystals as well (see e 



. g . Q,0 : 



17]). However, at frequencies where a cubic 



lattice becomes a photonic crystal the refraction index and the wave impedance both become 
strongly dependent on the propagation direction [22[]. This implies that effective material 
parameters e e s or/and /i fr of the crystal are very different for a wave propagating along 
one of the lattice axes and for a wave tilted with respect to this axis. This evidently means 
the absence of the isotropy. Of course, the random mixture of particles is more optically 
isotropic than the lattice with the same concentration of particles. However, in the present 
case the requirement of the negative magnetic response of the composite medium obviously 
implies the very dense packaging of magnetic scatterers. In the present case we cannot 
suppress the spatial dispersion using the randomization of the array. 

The comparison of the approximate dispersion diagram for the homogenized lattice with 
the exact dispersion diagram was done in some papers, e.g. in 2CJ? ]. In 20| the fee (face- 
centered cubic) and bec (body-centered cubic) lattices of metal-covered dielectric spheres 
were studied. It was found that in fee and bec lattices the spatial dispersion effects arise 
within the frequency band of the plasmonic resonance. The spatial dispersion in this region 
leads to closing the resonance band-gap with simultaneous strong anisotropy of the lattice 
refraction index. These effects correspond to the frequencies at which the distance between 
adjacent particles centers exceeds one quarter of the wavelength in the matrix. 

Unfortunately, in the known literature there are no general results for the frequency 
intervals where the effects of strong spatial dispersion arise in lattices. Of course, it is clear 
that the frequency range corresponding to D/X > 0.5 definitely implies the strong spatial 
dispersion since A = 2D is the Bragg resonance condition. In lattices of strongly scattering 
inclusions the condition of the absence of strong spatial dispersion becomes more restrictive 
since the polarization of inclusions shortens the effective wavelength A e fr compared to A. 



It is known that the artificial resonant magnetism in arrays of small complex-shape metal 
particles at microwaves practically corresponds to relation D p /X r > 0.1, where D p is the 
typical size of the particle 13j, |21[ and A r is the wavelength of the magnetic resonance of a 
particle in host medium. Notice, that the ratio D p /X ~ 0.1 . . . 0.15 is practically achievable, 
e.g. for well-known double split-ring resonators (SRRs). SRRs of this optical size can be 
considered as small enough to avoid the strong spatial dispersion in a compact lattice. Then 
we can consider the low-frequency stop-band of this lattice associated with the resonance of 
SRRs as the effect of the negative permeability. Really magneto-electric media are opaque 



when their material parameters have different signs / 



However, the metamaterial suggested in [14j corresponds to a quite large optical size of 
a magnetic scatterer at its resonance, where the condition D p /X r < 0.15 is far from being 
fulfilled. Really, consider the results for /i e g- illustrated by Fig. 2 (a,b) of [14J . The negative 
permeability was theoretically achieved at 630 or 680 THz with nanorings of the following 
geometry. Six or four metal nanospheres of radius a p = 16 nm are centered at the circle of 
radius Ro = 40 or i?o — 38 nm. This means that the size D p of the effective magnetic particle 
(i.e. the distance between the edges of opposite metal nanospheres) is 112 or 108 nm. The 
wavelength in the host medium with background permittivity = 2.2 at 630 THz is equal 
to X r = 323 nm. Then the size of the nanoring from 6 spheres is D p pa 0.35A r i.e. the optical 
size is rather large. The packaging of these nanorings in the host medium was assumed in [ijj] 
to be very dense, namely corresponding to their concentration N^r = (108)~ 3 nm -3 . This 
concentration is physically achievable for a bcc or a fee lattice of nanorings. However, even 
so dense packaging corresponds to the distance between the centers of adjacent nanorings 
D pa 0.4A r . Having in mind results of 2^] we believe that in this case the strong spatial 
dispersion arises and the isotropic magnetism is impossible. 

The clear way to decrease spatial dispersion in a lattice of magnetic scatterers i.e. to 
suppress the angular dependence of /x e ff is to decrease the size of a scatterer keeping its 
magnetic polarizability. In this way we will also decrease the optical distance between the 
centers of adjacent magnetic particles. Below we modify the design of magnetic particles 



suggested in 



14j in this way. We obtain the ratio D p /X r ~ 0.2 . . . 0.25 which is presumably 



below the range of the strong spatial dispersion. Simultaneously, we link our design to 
existing nanotechnologies. 
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II. CORE-SHELL MAGNETIC CLUSTERS 



The design of the optical magnetic scatterer suggested in this paper is based on the 
possibility to fabricate arrays of silver or gold colloidal nanoclusters located on spherical 
dielectric cores in a liquid host medium. The fabrication technology uses self-assembled 
adhesion of metal nanospheres on the silica core. The idea of such self-assembling aggregates 
was probably first realized in 



23]. However, this process led to the complete covering of 
the silica core by contacting gold nanoparticles. The process chemically stimulated by 
some acids which allowed one to control its speed and to obtain the needed number of 



nanocolloids per one core was described in 



241 ]. Later, nanoclusters comprising few tens of 



colloidal particles on a silica core of diameter of the order 100 nm were reported in works 
25)]. For obtaining the maximal electromagnetic response for given sizes (of the core and 
of the nanocolloids) their number N tot per one nanocluster should be as large as possible. 
However, it is difficult in this way to prevent the touching of nanocolloids which is destructive 
for their plasmonic resonance that will degrade due to the tunneling of electrons between 
nanocolloids. Practically, minimal separation d between nanocolloids should be comparable 
with the radius of colloidal particles a p . Below we respect the condition d ~ q „/2 that 
restricts N tot . This can be practically obtained using the technology described in [26( and 
2?]]. In these works the silica core of diameter 100 nm (the technology allows in principle 
to reduce this size to 20-30 nm) was covered by mutually touching polystyrene nanospheres. 
Instead of simple polystyrene nanospheres one can use core-shell particles with plasmonic 
nanocolloids as cores. This method will allow us to reliably control the separation d between 
colloidal particles which is simply equal to the double thickness of the polystyrene shell. 

Such a nanocluster is shown in Fig. [1] (a) (a general view) and (b) (a cross section). It 
should have the isotropic resonant response to the electromagnetic field of light. At one 
frequency the induced electric dipole dominates, and the polarization of colloidal particles is 
parallel to the electric field of light. At another frequency the induced magnetic dipole dom- 
inates, and the polarization of colloidal particles is azimuthal with respect to the magnetic 
field. In other words, the applied magnetic field form the effective nanorings of colloidal 
spheres glued to the silica core. Both electric and magnetic resonant frequencies origin 
from the plasmonic resonance of an individual colloidal nanosphere. In this paper we study 
only the magnetic resonance of such a nanocluster. Therefore we call it a magnetic nan- 





(a) (b) 

FIG. 1: A magnetic nanocluster (a bulk isotropic scatterer with the magnetic resonance): (a) - a 
general view of a cluster of metal nanocolloids with dielectric core, (b) - the resonant azimuthal 
electric polarization is excited in the discrete "shell" formed by plasmonic nanocolloids when 
the applied field is the optical magnetic field directed along z. 



ocluster (MNC). Bulk arrays of MNC can be created not only in liquid but also in porous 
matrices where the separation between the pores should be small enough to obtain the high 
concentration of MNC and a strong magnetic response of the medium. 

Below we will see that the frequency of the magnetic resonance of MNC is strongly 
reduced compared to the frequency Q p of the plasmonic resonance of a single nanocolloid. 
This effect results from the strong electrostatic coupling between the nanocolloids in the 



14[ , and it leads 



MNC. This electrostatic coupling exists also in a single nanoring studied in 
to a certain reduction of the magnetic resonance frequency compared to Q p . However, in 
the present geometry this reduction is much stronger since the most important electrostatic 
interaction is that between colloids having parallel electro-dipole moments. 

The total number N tot of colloidal nanospheres per one MNC can be controlled by the 
self-assembling technological process. Also, it can be adjusted by properly chosen initial 
concentrations of nanocolloids and silica spheres in the host liquid. N tot can be expressed 
through the core radius a, the colloid radius a p and the separation d between metal spheres 



(a) 



(b) 



FIG. 2: Two models of the magnetic polarization of MNC in the applied magnetic field directed 
along z: (a) - nanocolloids form regular rings around the z— axis, the colloidal particles are equidis- 
tant in the rings of radiuses Rj (induced dipole moments are shown by bold arrows) and the rings 
are equidistant; (b) - nanocolloids form around the core (shown as a dashed circle) an effective 
continuous shell with azimuthal bulk polarization Pa,. The shell of thickness 2a p + d can be char- 
acterized by resonant isotropic permittivity e^. 



Here [A] denotes the integer part of the number A. In this formula the curvature of the 
portion A of the spherical surface with radius Rq = a + a p + d/2 which is cut of the sphere of 
radius Ro (this sphere centers the nanocolloids and their dipole moments are located on it) 
by one core-shell particle of radius a p + d/2 is neglected. In other words, in our calculations 
the surface A of the effective sphere Rq per one colloidal particle is assumed to be a planar 
square with side 2a p + d. 

III. TWO MODELS OF MNC 

Two approximate models for obtaining the magnetic polarizability a mm of a MNC are 
illustrated by Fig. [2] (a) and (b). The agreement of their results can be hopefully considered 
as a validation of these results. These models describe in two different ways the electromag- 
netic interaction of colloidal particles. In the first model the nanocolloids on the dielectric 



(see Fig. HD: 




(1) 
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core are assumed to form regular rings around the z— axis. The rings of radiuses Rj are dis- 
tanced by a p /2 from one another. The distance between the nanoparticles in any ring is also 
equal to d — a p /2. The presence of the polystyrene shell of nanocolloids does not influence 
to the result, since the permittivity of the shell is equal to that of the matrix (e n = 2.2). 

In the second model the discrete structure of core-shell plasmonic nanoparticles on the 
central silica core is treated as a nanolayer of effective medium of thickness 5 = 2a p + d 
covering the silica core. The azimuthal electro-dipole polarization and the isotropic 
permittivity El of this layer can be easily calculated and we find a mm in the closed form 
calculating the magnetic moment of this homogenized shell. 

The presence of the central silica core also does not influence to the result due to the 
absence of its electric polarization by the applied magnetic field and smallness of its magnetic 
moment which is not resonant. It was checked analytically that the interaction of this 
small magnetic moment with electro-dipole moments of colloidal particles is negligible. The 
permittivity of the silica core e s = 4 > e n = 2.2 influences only weakly to the response of 
colloidal spheres. The polarizability a of a single colloidal particle is assumed to be the same 
as if nanocolloids were located in the uniform medium with permittivity e n '- 



a 



4irals e h - 



(2) 



"e m + 2e h J 6n n 
Here e m is the permittivity of the metal which is taken in the same form as in jl4j ]: 

F F (27T/ p ) 2 pKfUU 

<^m — fc i o r * 5 • \o ) 

U> 2 US 1 

For silver colloids we have the plasma linear frequency f p = 2175 THz, the damping fre- 



quency fa = 4.35 THz and t 



re ultraviolet permittivity = 5 



Q. 



As it was explained in [14J the magnetic moment of nanorings is excited by the magnetic 
field of the light wave. To calculate the magnetic polarizability of a MNC we assume that 
it is excited by magnetic field polarized along z [see Fig. [TJ (b)] that in the model of the 
magnetic excitation assume to be uniform across the circle of radius Rq and varying with 
time as exp(— tut)) with amplitude H Q . From Maxwell's equations this applied magnetic 
field H a is associated with the azimuthal electric field E a vanishing at the center of MNC: 

H a = z H , E a = <f (ir}H kR), (4) 

where rj = a/ Jio/so^h 1S the wave impedance of the host medium of permittivity s n - Unit 
vectors zq and </>o correspond to the cylindrical coordinates (R, <f), z) shown in Fig. [1] (b). In 



fact, the applied electric field given by (J4]) is associated with slightly non-uniform magnetic 
field H a = z H (l — k 2 R 2 /2), but this small correction to the applied magnetic field matters 
nothing for our model. We calculate the magnetic moment of MNC induced by H a whose 
amplitude at the center of MNC is H through the electric polarization of nanocolloids by 
the electric field associated with H a . 

The 0-oriented dipole moments pj of colloidal particles of any j-th ring are identical 
in polar coordinates. The z-directed magnetic moment rnj of the j-th ring can be found 
using formula (8) of [l4]. The total magnetic polarizability of MNC is then obtained as 



a m m = Y^ a mm, where a™L = mj/H . 

Notice, that the obvious (not sufficient) condition of the absence of spatial dispersion in a 
bulk array of MNC is the independence of the magnetic polarizability a mm of an individual 
MNC on the distribution of the applied magnetic field around the center of MNC, if the 
condition of the magnetic excitation is respected: H a (R = 0) = z H , E a (R = 0) = 0. This 
condition can be satisfied for example by a standing plane wave, i.e. a pair of plane waves 
travelling along the a;— axis: 

H a = z ^-(e-i kx + e***), E a = y^{e~ lkx - e jkx ). 

The comparison of a mm calculated for this azimuth-dependent excitation of MNC and that 
calculated for excitation (Tj0) would be a useful exercise, and we plan it for the next paper. 
It could indicate the frequency bounds of the spatial dispersion related to the finite size of 
MNC. However, in the present paper we restrict the analysis by a more simple excitation 
case @. 

In both models the reduction of the resonance frequency of MNC compared to the plas- 
monic resonance Q p of the individual colloidal particle in the host medium with Eh = 2.2 is 
determined by the electromagnetic coupling of colloids. In fact, this reduction will be even 
stronger than it is predicted below. This is because the presence of the silica core makes the 
effective host medium permittivity slightly larger than e^. We neglect this small effect for 
simplicity of the model, however, its role is positive. 



A. The model of regular rings 



Consider a dipole p = p(j)4>o shown in Fig. [T] (b) that belongs to the j-th effective ring 
of MNC. Assume that N tot ^> 1 (practically the model is applicable when N tot > 15 . . . 20). 
Then it is easy to show using auxiliary spherical coordinates and formula ([T]) for N tot that the 
radius of the j-th ring is approximately equal to Rj ~ Rq sin(2j ^/7r/N tot ). The z— coordinate 



of the j-th ring is z(j) ps R cos(2j^n/N tot ). The number of dipoles Nj in the ring is equal 
Nj 2TcRj/(2a p + d). The number of such rings in MNC is equal N r ps ^TiN tot /2. The 
cylindrical coordinates of the dipole p in the j-th ring are z = zj, R = Rj, = 2nq/Nj where 
q = . . . (Nj — 1) is the integer number determining the position of the dipole within the 
ring. 

The field produced by the dipole p at the point with coordinates (z f , R', 0') can be found 
from the standard formula: 



E(R, R') 



e ~jkr / e ikr ike ikr 

k 2 (r x p) x r — h (3rp ■ r — pr 2 ) I — — 



^3 \ i^*4 



(5) 



Ane e h 

Here r = ^/R 2 + R' 2 - 2RR' cos(0 - 0') + (z — z') 2 is the distance between radiating p(j) 
and receiving p(ri) dipoles of MNC. Due to the azimuthal symmetry of the problem only the 
0— th component of the field produced by all rings of dipoles p(j) is nonzero at the center 
of the receiving dipole p(n), i.e. at point (z', R', 0') . In other words the local electric field 
acting on the dipole p(n) is orthogonal to the vector R! shown in Fig. [T] (b). Respectively, 
only the azimuthal component of the vector E(R, R'), should be taken into account. 

From (jSj) the scalar interaction coefficient of dipoles p(j) and p(n) can be easily derived. 
It is defined as the 0-th component of the field produced by the unit azimuth-oriented 
dipole located at point (z = Zj,R — Rj, = 2irq/Nj) and calculated at point ( R! = 

i? n ,0' = 2vr S /A„): 

n ikr 



Q n l = 



(kr) 2 (R' 2 + (z- z')R' - RR' + cos(0 - <p')(R 2 + (z - z')R) - RR 1 ) 
3 J R J R'sin 2 (0-0')(l - ikr) -cos(0-0')(l - ikr)r 2 ] . (6) 



The dipole moment of the receiving dipole is equal to p — aE loc . Here the local field is the 
sum of the external electric field (0J E a (R') = (iH r]kR f /2) and all dipole fields: 



E loc = E a (R')+p(j)J2Q 



>> j 

qs ■ 



in 



This way we obtain the system of equations for dipole moments of colloidal nanospheres of 
any ring: 

9=1 j^n 9=0 

Here the term with Qq™ describing the interaction of the dipole p(n) located at (z 1 = z n , R' = 
R n , (J)' = 0) with other dipoles of the same n-th ring is shared out. The expressions for 
coefficients Qq™ entering this term and given by (J6j) with z = z' are described by (EH). There 
is no coefficient Qqq in (j7|) since the summation starts from q = 1 (see also in 141]). 

Solving the system (171) we find all dipole moments p(n). The magnetic moment of j-th 
ring is calculated as in [14 1: 

m . = -^(j)a^ (g) 

Then we obtain the magnetic polarizability of the MNC as the sum of m~ letting H = 1: 



j 

N r 

2 ~ 



limn, — — £ ;PU) N jRj- 
3=1 



The relative effective permeability of the composite medium is given by 14j : 



1 

^ = 1 + Api -i ~ - (9) 

ls MNC U mm 3 

Here Nmnc is the volume concentration of MNC that can be expressed though the effective 
volume per one magnetic scatterer N MNC = 1/Vo = 1/D 3 . In simple cubic lattices D > 
Dp = 2(a + 2a p ), and D is the unit cell size. However a more compact arrangement of MNC 
is also possible when D < D p and N MNC > V-^p- 

The inverse polarizability of a nanocolloid in (jTJ) corresponds to formula (j2J) and contains 
the term (— ik 3 /67reoEh) that describes the radiation damping. The radiation damping of the 
magnetic dipole with magnetic polarizability a mm should be described by the term (— z/c 3 /67r) 
[14J]. It is known that the radiation damping is cancelled out in regular 3D arrays and for 
regular lattices we should have used instead of (jUJ) the relation [l4|: 

^ = 1 + AT' 1 fa-^ + ^-l ' (10) 

^ MNC \ a mm ' 1 Q n ) 3 

However, the dissipative losses due to the plasmonic resonances of metal nanospheres 
strongly dominate over the contribution of radiation losses into the imaginary part of the 
permeability, i.e. the difference in results of (jUJ) and of ffTUl) is negligible. 

1 1 



B. The model of an effective shell 



Replacing the discrete plasmonic structure by the continuous shell we introduce the bulk 
polarization P^ that can be expressed through the averaged 4>— polarized electric field E &y 
in this shell and its unknown permittivity e^: 

P<p = e (e L - e h )E av (11) 

The field E av is related with the local field acting on any colloidal nanosphere by the Clausius- 
Mossotti relation: 

E av = E\ oc — — — . (12) 

Here V\ = 2a p (2a p + d) 2 is the volume per one colloidal nanosphere and p = P^Vi is the 
dipole moment of the reference nanosphere. Using the formula p = aE\ oc together with ifTTl) 
and ( fl2l) we come to the Lorentz-Lorenz formula: 

£ L = £ hi 1 + aeoel l j ■ ( 13 ) 

The definition of the magnetic moment of any volume V comprising polarization currents 
j reads as: 

1 

m = - / j x rdV. 
2 J J 

v 

It can be rewritten for the MNC in terms of the bulk polarization Pf. 



m = j P^RdV. (14) 

v L 

Here Vl = 4:nRl(2a p ) = 87ra p (a + a p ) 2 is the volume of the spherical layer with central radius 
Rq and the thickness 2a p . The integration of the bulk polarization across this layer can be 
replaced by simple product P^(2a p ). The polar radius R that enters (JT3J) can be expressed 
in spherical coordinates as R{9) = Rq sin 9. After substitution of ( fTTj) into ( fl4l) we obtain: 

2-7T 7T 

m = -icu£ a p (e L - e h )Rl J d<f) J d6 sin 9E av (9)R(9). (15) 

o o 

The averaged electric field of azimuthal polarization at the circle of radius R(9) is re- 
lated to the magnetic field H at the center of the MNC as E av (9) = iu){i R(9)H /2 = 
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(iHof]kRo sin 9/2), and after this substitution into (fl5|) and trivial integration we come (let- 
ting H — 1) to the following formula: 

a mm = m= —(k R ) 2 V (e L - e h ), (16) 

where it is denoted Vq = R^a p and k = k/ y/e~h = uo^/Eq^o is the free space wave number. 

Substituting ( f!3l) into (fl6l) we obtain the final closed-form formula for the magnetic 
polarizability of individual MNC: 

a mm = 47r(fco# )Vo - 3£o£ ^ (17) 
The permeability can be then found using ([91) or ( flOl) . 



IV. RESULTS AND DISCUSSION 

First, the results obtained in [jjj] were reproduced using the first model (the second model 
implies very large N tot and cannot be applied to the case of a single ring in MNC). Namely, 
the system ([7]) was solved for N r = 1, when it was taken also N tot = 4, Rq = 38 nm and 
a p = 16 nm. This geometry corresponds to the total size of the magnetic cluster D p = 108 
nm and to the spherical core radius (recall, that in the present theory the difference of the 
core permittivity from that of the host medium is not taken into account) a = Rq — a p = 22 
nm. The concentration Nmng of effective magnetic scatterers was assumed in [yj] to be 
equal Nmng — 95~ 3 nm -3 (almost touching MNC in a fee or a bec lattice), the metal of 
colloidal particles was silver with permittivity (J3]). 

The result for the effective permeability is presented in Fig. |3] (a). It reproduces Fig. 2 
(b) of [ijj] with high accuracy. This agreement can be considered as a check. 

If we keep the same separation of colloidal particles d = 10 nm as in the previous example 
the maximal possible number of colloidal particles for the spherical core of radius that 
corresponds to this example is equal to N tot = 16. It corresponds to N r = 3. Keeping the 
same concentration Nmnc — 95~ 3 nm~ 3 of same MNC we obtained with this geometry the 
result shown in Fig. [3](b). This result was also obtained using the first model (for N tot = 16 
the second model is inadequate). The reduction of the magnetic resonance frequency due 
to the presence of two additional nanorings compared to Fig. |3] (a) is not dramatic though 
visible. 
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(a) (b) 

FIG. 3: Effective permeability of the array of MNC with sizes D p = 108 nm, a = 22 nm, a p = 16 
nm (Ro = 38 nm) hosted in the matrix with Eh = 2.2. |(a)| One ring of nanocolloids in every MNC, 
the concentration of MNC is Nmnc = 95 -3 nm~ 3 . |(b)| Three nanorings in every MNC, other 
parameters are the same. Real and imaginary parts of the permeability are shown by solid and 
dashed lines, respectively. 



More promising results were obtained keeping the same total size of MNC D p = 108 nm 
for a p = 13 nm, d = 6 nm, that implies a = 22 nm. This geometry corresponds to one 
effective ring with 9 colloids, two rings with 5 colloids and two rings with 3 colloids. In this 
case the number of colloidal particles is large enough and the qualitative agreement of two 
models was obtained. The concentration of MNC in Fig. H] is the same as in the previous 
example: Nmnc — 95~ 3 nm' 3 . The two models show approximately the same resonant 
frequency and the same magnitude of the resonance. The second model ignores, of course, 
the higher magnetic resonances that should arise in interacting nanorings and are revealed 
if we use the fist model. However, these higher resonances are weak and have no practical 
meaning. 

The reduction of the resonant frequency in Fig. |4] compared to Fig. [3] means that the 
resonant size of the MNC reduces from 0.35A r as in the previous example to 0.25A in the 
present one. Here A r is the resonant wavelength in the host medium. 

High magnitude of the Lorentz resonance in Fig.HJmeans that we can also reduce the con- 
centration of MNC making the array of MNC easier for manufacturing. The results depicted 
in Fig. [5] (a) and (b) are obtained using the first model for Nmnc — 110~ 3 nm -3 instead of 
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108 nm, a = 22 nm, a p = 13 



FIG. 4: Effective permeability of the array of MNC with sizes D p 
nm, d = 6 nm hosted in the matrix with Eh = 2.2, the concentration of MNC is Nmnc = 95 -3 
~ 3 . (a) : First model of the interaction of nanocolloids. |(b)| Second model of the interaction of 



nm 



nanocolloids. 
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(a) (b) 
FIG. 5: Effective permeability of the array of MNC with sizes D p = 108 nm, a = 22 nm hosted 



in the matrix with = 2.2, the concentration of MNC is Nmnc = HO nm 3 . (a} 
d = 6 nm. |(b)[ a p = 12.5 nm, d = 6 nm. 



a p = 13 nm, 



Nmnc = 95~ 3 nm -3 as in the previous case. The difference between two plots in Fig. [5] (a) 
and (b) is determined by a 1 nm difference in the radiuses of plasmonic nanospheres, and 
demonstrates how sensible is the magnetic response of MNC to the deviation of parameters. 
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V. CONCLUSIONS 



In the present paper we developed the idea of the resonant optical magnetism in its 
isotropic variant modifying the known design of optical magnetic scatterers suggested in jjjl . 
The array of isotropic magnetic nanoclusters proposed in the present paper can be obtained 
in a liquid or porous matrix using the existing nanotechnologies. The metal nanocolloids 
are located in the core-shell particles attached to the silica core. Their electrostatic mutual 
coupling helps to reduce the resonant frequency for the given size of a nanocluster compared 
to the frequency of the plasmonic resonance of a single nanocolloid. This reduction is the 
essential result of the present paper since it allows one to decrease the spatial dispersion of 
the dense array (lattice) of nanoclusters. In other words, the suggested geometry is helpful 
for the isotropy of the effective resonant permeability of the array. The further optimization 
of the structure, further validation of the suggested model by full-wave numerical simulations 
are planned in the next paper. 
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